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In this work we study the well known contact term, which is the key element in resolving the 
so-called U{1)a problem in QCD. We study this term using the dual Holographic Description. 
We argue that in the dual picture the contact term is saturated by the D2 branes which can be 
interpreted as the tunnelling events in Minkowski space-time. We quote a number of direct lattice 
results supporting this identification. We also argue that the contact term receives a Casimir - 
like correction ~ (Aqcd-R)"^ rather than naively expected exp(— Aqcd-R) when the Minkowski 
space-time lZi,,\ is replaced by a large but finite manifold with a size R. Such a behaviour is 
consistent with other QFT-based computations when power like corrections are due to nontrivial 
properties of topological sectors of the theory. In holographic description such a behaviour is due 
to massless Ramond-Ramond (RR) field living in the bulk of multidimensional space when power 
like corrections is a natural outcome of massless RR field. In many respects the phenomenon 
is similar to the Aharonov -Casher effect when the "modular electric field" can penetrate into a 
bJ[) superconductor where the electric field is exponentially screened. The role of "modular operator" 

^ ■ from Aharonov -Casher effect is played by large gauge transformation operator T in 4d QCD, 

resulting the transparency of the system to topologically nontrivial pure gauge configurations. We 
discuss some profound consequences of our findings. In particular, we speculate that a slow variation 
of the contact term in expanding universe might be the main source of the observed Dark Energy. 
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I. INTRODUCTION 



The holographic picture of QCD is known to provide some deep insights into the strongly coupled dynamics. 
The most well known example is the dual model of gluodynamics at zero and nonzero temperatures [1], and its 
generalization [2] when light fermions are included into the system. 

The goal of thsi work is to study the well known contact term which is the key clement in in resolving the so-called 
U{1)a problem in QCD [3-5], see also [6-8]. We want to understand some deep properties of the contact term from 
the dual perspective. 

The basic tools in this study will be the D branes which are the crucial elements of the dual description. It is 
well known that DO brane extended along X4 can be identified with conventional instanton [9]. At the same, the D2 
brane wrapped around both periodic coordinates was identified with magnetic string [10], while D2 brane wrapped 
around X4 was identified with D2 domain wall [11]. These two branes will play an important role in our studies. For 
completeness, we should also mention the configuration of the D6 brane wrapped around the compact 5^ part of the 
dual geometry was interpreted as the domain wall separating two vacua in pure gluodynamics [12]. The corresponding 
generalization with quarks was given in [13, 14]. A large number of other D defects from this Zoo have been discussed 
in [15]. 

The crucial element which helps to make the identifications is the study of the confinement- dcconfincment phase 
transition which is interpreted in the dual picture as the Hawking-Page phase transition [1]. In the dual picture the 
wrapping around X4 is stable above the phase transition T > Tc a,s a, result of cylinder geometry, while it is unstable 
below the critical temperature T < Tc as a result of cigar type geometry, see definitions and details below. The 
wrapping around r does the opposite: the corresponding D brane is stable at small temperatures T < Tc and unstable 
at high temperatures, T > Tc- Since the contact term is sensitive to the ^-parameter we shall also discuss the 9 
dependence from the dual perspective. 

To formulate the problem we start from the conventional definition of the topological susceptibility x which plays 
a crucial role in resolution of the U{1) problem [3-8] 



lim / d4xe'^-(T{g(a:),g(0)}), (1) 



where 9 is the 9 parameter which enters the Lagrangian along with topological density operator q{x), see precise 
definitions below. 

It is important that the topological susceptibility x does not vanish in spite of the fact that q — df^K'^ is total 
divergence. Furthermore, any physical asymptotic states gives a negative contribution^ to this diagonal correlation 
function 

Xdispe.ive - lim / d4xe^'^-(T{g(x),g(0)}) ^ lim ^"''^'f ^^'f^ -M! < Q. (2) 

where mc is the mass of asymptotic state, fc — >■ is its momentum, and (OlglG) = cg is its coupling to topological 
density operator q{x). At the same time the resolution of the U {1)a problem requires a positive sign for the topological 
susceptibility (3), see the original reference [5] for a thorough discussion, 

Xnon-di.pcrsivc = Hm / x^'^^T {q{x) , g(0)}) > . (3) 

fc— >0 J 

Therefore, there must be a contact contribution to x, which is not related to any propagating physical degrees of 
freedom, and it must have the "wrong" sign. The "wrong sign" in this paper implies a sign which is opposite to any 
contributions related to the physical propagating degrees of freedom. In the framework [3] the contact term with the 
"wrong sign" has been simply postulated, while in refs.[4, 5] the Veneziano ghost had been introduced to saturate the 
required property (3). These two descriptions are equivalent as they describe the same physics. 

It interesting to note that the "wrong" sign in topological susceptibility (3) is not the only manifestation of this 
"weird" unphysical degree of freedom. In particular, one can argue that this term also contributes with a wrong sign 
into the entropy. While the entropy itself is a positively defined entity, the corresponding gauge invariant contribution 
related to the topological susceptibility contribute with the negative sign [16, 17]. 

The goal of this paper is to investigate this "weird term" using the holographic description. We also want to study 
the behaviour of the system when size of then system is large but not infinitely large. The last property, as we 



^ We use the Euclidean notations where path integral computations are normally performed. 
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argue below, might be relevant for cosmological applications, see section V. We also want to study the correlation 
between drastic changes in 9 behaviour when the phase transition is crossed at T^- This point has been emphasized in 
[9, 10, 15, 18] from holographic perspective as well as from quantum field theory (QFT) viewpoint. This correlation 
is also supported by direct lattice studies, see review [19] and references on the original papers therein. 

The paper is organized as follows. In section II we overview the nature of the contact term from QFT viewpoint. 
In Section III we describe our model based on the Nc D4 branes with one compact world-volume coordinate. We also 
overview some D- defects relevant for our analysis. We treat D defects in the probe approximation when they do not 
deform the dual geometry. In Section IV we formulate our proposal on contact term from holographic perspective. 
We present a number of arguments supporting our construction, including the comparison with the direct lattice 
measurements. We also argue in section IV B that the contact term receives a Casimir -like correction ~ (Aqcd-R)~^ 
rather than naively expected exp(— Aqcoi?) when the Minkowski space-time TZ^^i is replaced by a large but finite 
manifold with a size R. In holographic description such a behaviour is due to massless Ramond-Ramond (RR) field 
living in the bulk of multidimensional space when power like corrections is a natural outcome of massless RR field. 
In many respects the phenomenon is similar to the Aharonov -Casher effect when the "modular electric field" can 
penetrate into a superconductor where the electric field is exponentially screened. In Section V we outline some 
profound consequences of our findings. In particular, we speculate that a slow variation of the contact term in 
expanding universe might be the main source of the so called Dark Energy (DE). 



II. CONTACT TERM IN QCD FROM QFT VIEWPOINT. 



In this section we want to provide some intuition on the nature of the contact term using some model computations. 
The simplest possible example is two dimensional QED2 in the Kogut-Susskind formulation [20] where all essential 
elements related to the contact term can be explicitly worked out as discussed in refs. [16, 21]. However, for purposes 
of the present work we present the relevant exact results for four dimensional "deformed QCD" formulated in [22] 
where all computations can be explicitly performed as the model is in a weak coupling regime, see next subsection II A. 
We reformulate the same results for strongly coupled 4d QCD using effective description in terms of the Veneziano 
ghost as originally developed in [4, 5], see section II B. 



A. Contact term in weakly coupled "deformed QCD" 



In the deformed theory an extra term is put into the Lagrangian in order to prevent the center symmetry 
breaking[22]. It can be arranged in such a way the constructed theory remains confined at high temperature in 
a weak coupling regime, such that there is no order parameter to differentiate the low temperature confined regime 
from the high temperature confined regime. Therefore, one can use this weakly coupled gauge theory to test some 
deep theoretical ideas about confined phase. 

The topological susceptibility in this model can be explicitly computed and is given by [23] 



Xqcd 



d^x{q{^)q{0)) 



C 



NrL 



(5(x) -m^,- 



(4) 



where the monopole fugacity C can be explicitly computed in this model [22]. The monopoles are pseudoparticles 
in this model. Therefore, the monopole fugacity C should be understood as number of tunnelling events per unit 
time per unit volume. There are two terms in formula (4). The first term represents the non-dispersive contribution 
not related to any physical propagating degrees of freedom, while the second terms represents the conventional 77' 
contribution. As one can see the first term contributes with sign plus while the second term contributes with sign 
minus in accordance with general equation (2). Another important lesson from this equation: the Ward Identity (WI) 
expressed as XQCoijUq = 0) = is automatically satisfied as a result of cancellation between the two terms. 

It is important to note that the number of tunnelling events per unit time per unit volume (4) in pure gauge theory 
in this model (with no quarks) is related with the absolute value of the energy density of the system. Indeed, 



NcC 



■ cos 



Nr. 



c 



(5) 



where we keep only the lowest branch Z = in expression for cos ^^^^ j to simplify formula (5), see detail discussions 
with complete set of references on this matter in [23]. In different words, the contact term in pure gauge theory 
Xym = ''^^ ^® interpreted in terms of number tunnelling events between different topological sectors in the 
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system. Therefore, there is no surprise that it has "wrong sign" as the relevant physics can not be described in terms 
of propagating physical degrees of freedom, but rather, is described in terms of the tunnelling events between different 
(but physically equivalent) topological sectors in the system. 

B. Contact term in terms of the Veneziano ghost in strongly coupled QCD 

Our discussions thus far were limited to the model where all computations are justified as a result of the weak 
coupling regime enforced by deformation [22]. Now we want to explain the same physics, the same contact term with 
"wrong sign" but in different way. We want to formulate this physics in terms of the Veneziano ghost which effectively 
describes the dynamics of the topological sectors in strongly coupled regime. Such an alternative description in terms 
of the Veneziano ghost [4, 5] will play an important role in our identification of the relevant configurations in next 
section when we will discuss the contact term using the dual holographic description. 

The topological susceptibility in the chiral limit can be easily computed in the model [4, 5] using the Veneziano 
ghost and it is given by [24, 25]: 

d^x{T{q{x), g(0)}> = • J d^x [5\x) - ml,D^{m,,x)] . (6) 

The Green's function D'^^mri'x) in this expression describes free (in the chiral limit) massive 77' field and satisfies 
standard normalization J d'^xmfy D'^{mjjix) = 1. 

The structure of this expression (6) is identical to our formula (4) computed in weakly coupled "deformed QCD" . 
In particular, the term proportional —D'^{mr,'x) with negative sign in eq. (6) is resulted from the lightest physical 77' 
state of mass to,,'. At the same time, the contribution with "wrong sign" expressed by 5'^{x) in eq. (6) is nothing else 
but the contact term 

r f^,TO^, r 

XYM ^ j d^x{T{q{x), qm)YM = • j d^x [5\x)\ . (7) 

In the model [4, 5] this term is saturated by the Veneziano ghost contribution. The required "wrong sign" for this 
contribution is due to the negative sign for the kinetic term for the Veneziano ghost. In our weakly coupled "deformed 
QCD" when all tunnelling processes can be explicitly computed this contact term in eq. (4) is expressed in terms of 
the monopole's fugacity NcC/L, and describes tunnelling transition between the topological sectors of the theory. It is 

natural to expect that the same interpretation should be also applied to dimensional factor ^' ^ which enters (7). In 
different words, this factor should be also interpreted in terms of tunnelling of some "objects" . However, the nature 
of these objects can not be easily identified in contrast with weakly coupled "deformed QCD" where the relevant 
pseudo-particles are obviously the monopoles saturating the contact term. We shall see however in next section that 
such an interpretation still can be advocated even in strongly coupled gauge theory. However, our arguments will be 
based on dynamics of "dual objects" rather than on dynamics of the original gauge fields. 

One should say that the main features represented by eqs. (6), (7) such as singular behaviour of the contact 
term with positive sign, and a smooth behaviour of conventional dispersive contribution are explicitly seen in lattice 
simulations [26-29]. A similar structure has been also observed in QCD by different groups [30-34] and also in 
two dimensional CP^~^ model [35]. Most important part for us from these numerical studies is that the singular 
behaviour of the contact term is not an artifact of any approximation, but an inherent feature of underlying gauge 
theory. Furthermore, there are no any physical scale factors (such as Aqcd) which would determine the singular 
behaviour of this term. In different words, the Veneziano ghost does model the crucial property of the topological 
susceptibility related to summation over topological sectors in gauge theories. This feature can not be accommodated 
by any physical asymptotic states as it is related to non-dispersive contribution with "wrong sign" . 

This concludes our discussions of the contact term within QFT approach. The corresponding computations give 
us a hint on the nature of the contact term, and quantum configurations which saturate it. These insights will play 
an important role in our analysis of the contact term in holographic description, where we attempt to identify the 
relevant configurations which saturate xym in the dual picture. 



III. DESCRIPTION OF THE MODEL 

We start with the Witten's background [1] for SU{Nc) Yang-Mills theory in large N,. limit, which is the gravity 
dual to the theory on D4 brancs wrapped around a compact dimension S^. In the supergravity approximation 
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the geometry looks as Miq = i?3.i x D x and the corresponding metric in notations [9] reads as 

ds' = i^f/^-dt' + 6,,dx^dx^ + f{u)dxl) + (^)-3/2(^ + u'dQl) 
Ho Ho ,t(u) 



m = l-i^f, Ro = {.g.N^f'\ R=^^{^f'\ (8) 
Hq u o ua 



The couphng constant of Yang- Mills theory and the radius of the compact dimension R are related as follows 

In formula (8) the confined phase corresponds to the geometry of a cigar in the (11^x4) coordinates with the tip at 
u = u\. At the same time, in de-confined phase the (u, X4) coordinates exhibit the cylinder geometry. The situation is 
reversed when one considers (r, u) plane instead of (x4, u) plane, where r is the Wick-rotated time coordinate r = it, 
T (X T + (3. It has been argued in [1] that the phase transition occurs exactly when one geometry replaces another. 
Important feature of this transition is as follows. The wrapping around the internal X4 circle is topologically stable, 
while the wrapping around the Euclidean time coordinate t is unstable above the phase transition point T > Tc- The 
opposite patter occurs of the two wrappings below the phase transition at T < Tc- Namely, the wrapping around 
the internal X4 circle is topologically unstable, while the wrapping around the Euclidean time coordinate r becomes 
stable at T < Tc. 

We want to discuss the topological features of the theory related to the 6'-dependence. The corresponding bulk field 
that couples to the topological density operator is the Ramond-Ramond (RR) field. The precise relation between 9 
parameter and RR field is known, 

9 = y J dx4 Ci (10) 

where Ci is the RR 1-form field. Important consequence of the relation (10) is that the 9 dependence of the world- 
volume theories on the D branes correlates with the wrapping around X4 coordinate. In different words, D brane 
configurations might be relevant for study the 9 dependence if the corresponding construction includes the wrapping 
around X4 coordinate. Otherwise, the D brane configurations can not be sensitive to the 9 parameter. 

We conclude this section with short description of D branes which will play an important role in our discussions 
which follow, see [15] with description of other members of this Zoo: 

DO instantons . 

The DO brane extended along X4 was identified as the YM instanton [9]. At T > Tc this wrapping is topologically 
stable, while at T < Tc the DO brane tends to shrink to the tip where its tension vanishes. Nevertheless, it has been 
argued in [15, 18] that the DO branes might be still important configurations at T < Tc- 
D2 string. 

The magnetic string is the probe D2 brane wrapped around Si parameterized by X4 and its tension is therefore 
proportional to the effective radius R{u) [10]. At small temperatures this wrapping is topologically unstable and the 
D2 brane tends to shrink to the tip where its tension vanishes. The magnetic string carries nontrivial 4d topological 
charge. 

D2 domain wall . 

The D2 brane wrapped around X4 was identified as D2 domain wall [11]. At small temperatures this wrapping is 
topologically unstable and the D2 branes tend to shrink to the tip where their tension vanish. At T > Tc the D2 
branes arc topologically stable objects and have been identified with Zjsi domain walls which are present in the system 
in high temperature phase. 
D6 domain wall . 

If D6 brane wraps around t it behaves as the domain wall which is a source of the corresponding RR-form. Such a 
configuration has been interpreted in ref. [12] as the domain wall which separates different metastable vacua known 
to exist in gluodynamics at large Nc- Its worldvolume theory on the domain wall has no 9 dependence as D6 brane 
wraps around r rather than around X4 which was the case for DO, D2 branes mentioned above. 



IV. HOLOGRAPHIC DESCRIPTION OF THE CONTACT TERM 



The computation of the contact term in holographic picture is known [12, 36], see also review paper [37]. The 
contact term has a positive sign and has expected iV" behaviour in accordance with QFT expression (7). The final 
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formula in holographic computations is expressed as a surface integral in 5-th dimension and in notations of section 
III is given by [12, 36, 37]: 

poo 

XYM - / dudu ^ M^A^iVO, h{u) ^ f{u) = 

J UA 

where M = dx4]~'^ is the only dimensional parameter of the problem, to be identified with Aqcd. Function h(u) 
in (11) is 4 — th component of the R-R field Ci which satisfies classical equation of motion 

[V5.9"^"^5'"'a./i(zi)] = 0. (12) 

The solution of this equation is given by expression (11), see [12, 36, 37] for the technical details. 

We have nothing new to add to those computations. The goal of this section is in fact quite different: we want to 
understand the physics which is hidden behind these computations. In different words, we want to understand the 
relation between the holographic formula (11) and the QFT based computations (6, 7). As we mentioned previously, 
the QFT based computations (6, 7) are supported by the lattice studies [28, 30-32], including the singular behaviour 
of the contact term with vanishing size in the continuum limit. This singular behaviour, as it is known, leads to a 
finite integral contribution to xym with the "wrong sign" in agreement with eqs. (6, 7). The understanding of this 
singularity in the dual description will be the crucial element in our analysis on a possible variations of the contact 
term when the background slightly varies, and when formulae similar to (11, 12) are not presently available. Our 
findings regarding these possible variations with slight variation of a background will play a key role in the application 
considered in section V where we interpret a tiny deviation of the 9— dependent portion of the vacuum energy (coded 
by Xym) in expanding universe as a main source of the observed dark energy. 



A. Emergence of the contact term as a tunnelling process. Percolation of D2 branes. 

The presence of non-dispersive contact term (7) in topological susceptibility in pure gluodynamics obviously implies 
that there must be unphysical pole in gauge variant correlation function 

X^YMik) ~ d^xe^''^TK^{x),K''m ^ xym - k^K ■ X^m (13) 

While the correlation function Xya/(^) is gauge variant itself, the merely presence of the pole, its position, its residue 
with "wrong sign" are physically observable parameters as corresponding correlation function (1), (7) is gauge invariant 
object. What is the nature of this pole? What it could be? How it can be interpreted from the holographic perspective? 
Normally, a pole at zero mass corresponds to a massless gauge boson. Or it might be a result of spontaneous symmetry 
breaking effect. What is a symmetry which could be responsible for behaviour (13)? Furthermore, as explained in 
section II this pole must have a residue with a "wrong sign" such that it can not be identified with any physical 
propagating massless degree of freedom. In QFT description this pole corresponds to the contact term which emerges 
due to the presence of different topological sectors in the system as argued in section II. In particular, in weakly 
coupled "deformed QCD" the contact term is saturated by monopoles which describe the tunnelling between those 
topological sectors as discussed in II A. In strongly coupled regime analytical computations based on underlying gauge 
fields are not reliable. However, one can argue that the contact term is saturated by the Veneziano ghost as reviewed 
in section II B. In this framework the Veneziano ghost should be treated as an auxiliary unphysical degree of freedom 
which does not belong to the physical Hilbert space. Nevertheless, it saturates topological susceptibility (6) and 
describes the dynamics of these degenerate topological sectors of the theory. 

The QFT- based insights presented above suggest that the origin for the contact term in holographic description 
should be also interpreted as a result of tunnelling events rather than presence of any massless propagating degrees of 
freedom. The tunnelling events in path integral approach in four dimensional Euclidean space should be manifested as 
a presence of some "objects" which live in four dimensions. We argue below that we can identify these objects as D2 
branes wrapped around 5*1 parametrized by Xi. As we mentioned in section III the tension for the D2 brane vanishes 
as a result of cigar geometry such that the D2 brane configuration shrinks to the tip. Therefore, an arbitrary large 
number of D2 branes could be produced, resulting in their condensation. The term "condensation" we use here is 
actually a jargon, as D2 branes are 4d time- dependent objects rather than 3d static objects when term "condensation" 
is normally used. The "tunnelling" or "4d -percolation" probably are more appropriate terms to describe the relevant 
physics in this case. However, we keep the term "condensate" as it is an appropriate term when the system is viewed 
from 5d perspective, which is precisely what holographic description is about. We shall present below a number of 
arguments supporting this identification. 

To begin with our arguments, we remind the reader the Bloch's theorem in condensed matter physics which states 
the following. Consider a system (such as a perfect lattice) which is invariant under displacement x ^ x-\-a symmetry. 



7 



The system in such circumstances can be described by the bands of allowed and forbidden energies classified by a 
quasi-momentum. Furthermore, a quasi-electron in allowed band is characterized by ungapped dispersion relation 

e(fc) = 2^7^ such that the corresponding correlation function has a pole at e(fc) — = while all information 

about tunnelling properties are coded by parameter m* and by residue of the Green's function. The presence of 
such pole implies that the corresponding quasi-particle can freely propagate to arbitrary large distances without any 
disturbance from the strongly interacting lattice as long as the symmetry x ^ x + a is respected and no any lattice's 
defects are included into the consideration. It should be contrasted with a case when an electron is placed into a 
system with a conventional potential barrier, in which case the correlation function will be characterized by a gapped 
dispersion relation with a pole at [e(fc) = + such that a particle can propagate only for a short time ^ uj^ . 

Our system with winding states quasi-momentum and large gauge transformation operator T with property 
T\n) = \n + 1) is analogous to the Bloch's displacement symmetry x ^ x + a with the "only" difference that the 
states \n) are not physically distinct states, and there is no any real "physical" degeneracy in the system. The physical 
vacuum state in QCD is unique and constructed as a superposition of \n) states. Nevertheless, the presence of the pole 
in eq. (13) can be interpreted as a result of tunnelling between topological sectors similar to the tunnelling events 
in the Bloch's system. It also explains the "wrong" sign in residues of the correlation function (13) as we describe the 
tunnelling of the Euclidean D2 objects between winding states \n) rather than the tunnelling of conventional physical 
degree of freedom between distinct vacuum states in condensed matter physics in the Bloch's case, when the residue in 
the Green's function (coefficient in front of the pole) is precisely the probability of the tunnelling. It further supports 
our interpretation of the coefficient /^/™^//4 in eq. (7) in strongly coupled QCD as a probability of tunnelling of 
"some objects" which we identify now with D2 branes. 

Now we can answer the question formulated in the beginning of this section: what is a symmetry which could be 
responsible for behaviour (13)? It is the invariance under large gauge transformations defined by the operator T 
such that [T, H] = 0, see Appendix A for a technical details. We come back to this interpretation at the end of this 
section when we discuss analogy with Aharonov -Casher effect and similarity between the large gauge transformation 
operator T and their "modular operator" commuting with the hamiltonian. 

The tunnel-based interpretation unambiguously implies that the corresponding D-defect which could be potentially 
responsible for tunnelling must be constructed with wrapping around Si parametrized by x^ rather than by r coor- 
dinate. In this case the corresponding D defect becomes the 6- dependent object (and therefore becomes sensitive to 
the transitions \n) |n + 1)). This wrapping around 5*1 parametrized by X4 leads to vanishing action at the tip of 
the cigar geometry at low temperature. The D2 branes responsible for the transitions |n) — >■ \n + 1) precisely satisfy 
both these properties, which explains our choice for a relevant configuration in the dual description. As we mentioned 
at the end of section III one type of the D2 branes was identified earlier in ref. [10] as magnetic strings. Another 
type of the D2 branes was identified as D2 domain walls [11]. The D2 magnetic strings carry 4d topological charge by 
construction, while D2 domain walls carry Zn charge [11]. Nevertheless, these D2 domain walls in combination with 
DO branes in strong coupling regime at T < Tc are capable to carry 4d topological charge. Indeed, the corresponding 
D0-D2 configurations can be explicitly constructed in some supersymmetric models [38] where analytical QFT based 
computations can be also performed. Similar D0-D2 configurations may also emerge in strongly coupled QCD at 
T < Tc as argued in [15]. 

Therefore, in what follows we shall not discriminate between these two options in our discussions as lattice studies 
can not presently precisely determine the dimensionality of fluctuating quantum objects. In particular, in ref. [34] it 
has been argued that these long range configurations actually might be characterized by Hausdorff dimension which 
gradually varies with cooling procedure. Therefore, we shall use a single term "D2 brane" to describe these objects 
in spite of the fact that these configurations are physically different. We expect that some future, more refine lattice 
studies will allow to discriminate between those two options (or rule out both of them)^. 

The recent Monte Carlo studies of Yang Mills theory have revealed a laminar structure in the vacuum consisting 
of extended, thin, coherent, locally low-dimensional sheets of topological charge embedded in 4d space, with opposite 
sign sheets interleaved, see original QCD lattice results [26-29]. One should add that a similar structure has been 
also observed in QCD by a different groups [30-34]. It has been also observed in two dimensional CP^~^ model [35]. 
We do not make any attempts in this work to cover this subject which has a number of subtle points. Instead, we 
concentrate on very few key properties of the gauge configurations which apparently make crucial contributions to 
the 6 dependent portion of the energy expressed in terms of the topological susceptibility (7). One can argue that 
such a complex structure observed on the lattice can be thought as a result of dynamics of the D2 branes described 
above'^ . 



^ I am thankful to Ivan Horvath for his comments on present measurements of dimensionality of the objects observed on the lattice 

simulations, and on future progress which should be expected as a result of these measurements. 
^ Similar idea with interpretation of the D6 branes in holographic description being responsible for topological susceptibility was re- 
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Indeed, it has been observed that the tension of the "low dimensional objects" vanishes below the critical temper- 
ature and these objects percolate through the vacuum, forming a kind of a vacuum condensate. Furthermore, these 
"objects" do not percolate through the whole 4d volume, but rather, lie on low dimensional surfaces 1 < c? < 4 [27]. 
The total area of these surfaces is dominated by a single percolating cluster of "low dimensional object". From 
holographic perspective these Monte Carlo results fits very nicely with conjecture that the observed structure can be 
identified with the D2 branes described above. In particular the tension of the D2 branes vanishes in confined phase 
as a result of cigar geometry in the holographic 5d description such that the D2 brane shrinks to the tip as we already 
mentioned. Therefore, an arbitrary large number of D2 branes could be produced, resulting in their condensation. 
The observed percolation (condensation) of these objects in lattice simulations unambiguously imply that they must 
have vanishing effective tension. Otherwise, only finite, not infinite percolating clusters could be observed, in contrast 
with observations. 

Another measurement suggests that the contribution of the percolating objects to the topological susceptibility xym 
has the same sign compared to its total value. Furthermore, both contributions have opposite sign in comparison 
with any propagating physical states (2), (6). These lattice studies are consistent with non-dispersive nature of D2 
branes, which can not be related to any propagating physical degrees of freedom. It unambiguously implies that the 
D2 branes should be thought as configurations describing the tunnelling events in Minkowski space, rather than as a 
real physical domain walls. 

Furthermore, the width of the "objects" apparently vanishes in the continuum limit. These measurements, again, 
support the tunnelling interpretation of the contact term. Indeed, it is well known property of conventional quantum 
mechanics that the individual photons penetrate an optical tunnel barrier with an effective group velocity considerably 
greater than the vacuum speed of light [42], see also review [43]. In different words, the tunnelling time can not be 
distinguished (experimentally) from zero. In quantum field theory context it means that the configurations responsible 
for tunnelling events must have vanishing size (in continuum), which is precisely what has been measured on the 
lattices. The same feature can be also seen on Fig 5 of ref.[32] where the contact term with vanishing size (in 
continuum) has a "wrong sign" in Euclidean lattice simulations. This again supports the tunnelling interpretation of 
this term in Minkowski space. 

B. Variation of the contact term with small variation of the background 

Our previous section IV A was devoted to the discussions of the contact term in QCD from holographic perspective. 
Essentially, we have not produced any new results in that section as the presence of the contact term is well established 
phenomena (supported by the lattice studies). We simply suggested a new interpretation of this term. However, the 
new interpretation may lead to some new physical results as we shall argue below. 

Main question we address here is related to a possible variation of the contact term when the background is slightly 
changed. For example, we wish to know what happens when infinite time coordinate r e ("00,-|-cx)) is replaced by 
a finite size ring r S (0, 27ri?), or what happens when the Minkowski space-time i?3^i is replaced by FRW metric 
characterized by dimensional parameter R ^ H^^ describing the size of the visible universe with H being the Hubble 
constant. 

Normally, one should not expect any strong dependence on very large distances as QCD has a mass gap ~ Aqcd- 
Therefore, one should expect the exponential suppression ^ exp(— Aqcd^) for any local observables. Formally, this 
expectation follows from the dispersion relation similar to (2) written in coordinate space which explicitly shows 
an exponentially weak sensitivity ~ exp(— Aqcd^) to large distances. The main point of this paper is that along 
with conventional dispersive contribution (2) there is also non-dispersive contribution (3) which is not related to 
any physical propagating degrees of freedom, and which can not be computed using the dispersion relations. This 
contribution as we discussed above emerges as a result of topologically nontrivial sectors in four dimensional QCD, 
and it may lead to a power like corrections with power p> \ rather than exponential like exp(— Aqcd-R)- In fact, 
this term in our framework is described by topologically protected massless ghost field (13) as discussed in section 
IV A. Therefore, in this framework, it is indeed quite natural to expect a power like corrections R~'p. 

We want to get some insights on this problem from the holographic perspective. We start our discussions with few 
general remarks. In principle, when we slightly change the geometry, one should compute the corresponding changes 
in the vacuum energy as a result of changes of the boundary conditions. The corresponding changes represent the 



cently advocated in rcf. [39]. One should remark that these D6 domain walls are in fact 6-independcnt configurations according to the 
classification presented in section III, and they can not interpolate between different topological sectors |n) and |n-|- 1). In different 
words, these D6 domain waUs interpolate between the ground state and an excited local minimum characterized by a small extra energy 
density AE ~ 1/Nc- These domain walls are similar to the physical static domain walls in ferromagnetic systems interpolating between 
physically distinct vacua. These D6 domain walls have finite action and finite width in contrast with thin tensionless D2 branes. Such 
domain walls in QCD within QFT approach in the presence of light quarks in large Nc limit when the -q' is light m? , ~ 1/Nc had been 
previously studied in [40]. 
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Casimir vacuum energy in multidimensional space due to new specific boundary conditions which were trivial ones 
in previous computations (11), (12) with infinite Minkowski space-time background. Technically, it could be very 
challenging problem as it requires the computation of a corresponding multidimensional Green's function in curved 
background which satisfies specific boundary conditions. However, on a general ground we expect that the contact 
term (11) receives a power-like correction ^ e^, not exponentially small corrections ^ exp(— 1/e) if deformation of the 
geometry parametrized by small parameter e ^ (_RAqcd)~^ ^ 1- 

Indeed, the computations of the Casimir vacuum energy in 5-dimensional space {x^^^u) can be accomplished by 
using effective Lagrangian describing the massless RR (10) or axion field"* living in the bulk of multidimensional space. 
This field Ci is responsible for the contact term (11). More importantly, the axion field is massless in multidimensional 
space as a consequence of corresponding gauge invariance. Therefore, while a computation of the Casimir vacuum 
energy in 5-dimensional space {x^,u) could be very challenging technical problem, an estimation of this energy is 
simple: we expect a power -like correction to the vacuum energy as a natural consequence of the dynamics of a 
massless axion field living in the bulk of multidimensional space, which is responsible for the contact term. 

We will test this argument (on power like correction) below by replacing infinite Euclidean time coordinate r S 
(— cxj,-|-oo) by a finite size ring t € {0,2ttR). As is known, sufficiently small parameter R corresponds to large 
temperatures T ^ where QFT based computations arc justified. After that we apply the same argument in strong 
coupling regime where QFT based analytical computations are not available. Nevertheless, we will quote some recent 
lattice numerical results supporting our argument. 

First, we consider de-confined phase when the length of the ring S]. is quite small and order of T^^. The contact 
term is expected to be exponentially suppressed at T > Tc in de-confined regime as DO brane wrapping around X4 
has a finite tension, saturates the topological susceptibility, and leads to the exponentially small contact term [9, 18]. 

This picture can also be easily understood from QFT viewpoint [15, 18]. Indeed, the wrapping around X4 corre- 
sponds to the well defined small instanton and one can use the standard instanton calculus to estimate the critical 
temperature Tc and the 9 behaviour above Tc'. 



Vinsti9,T >Tc)^ cos 9- e- 



1 > 



Xym{T > Tc) 



d^Vinst{9,T) 
89^ 



T-Tc 
Tc 



0, 



» 1/Nc 



1, Nc » 1. 



(14) 



Such a behaviour implies that the dilute instanton gas approximation at large Nc is justified even in close vicinity of Tc 
as long as "^j"^" :$> Such a sharp drop of the topological susceptibility xvAiiT) is supported by numerical lattice 
results [19] which unambiguously suggest that the topological fluctuations are strongly suppressed in de-conflned 
phase, and this suppression becomes more severe with increasing Nc- The behaviour (14) in de-confined phase should 
be contrasted with formula (11) corresponding to the confined phase when xym{T = 0) A^". 

Now we want to address the following question: how docs the contact term xy m vary when the radius of the ring 
slightly changes? We want to study these changes by remaining in a deep de-confined phase (T ^ Tc), such that 
no drastic variations are expected to occur. One can easily answer this question by noticing that small variation in 
T — > T + At corresponds to a small variation of the temperature T ^ T — AT with AT = (27rAT)~^ while keeping 
Tc unchanged. Therefore, a small variation of the background r — >■ r At results in a small (power-like) variation of 
the contact term 



Xym{T-AT) 



-aNa( — 



■-) 



Xym{T) 



1 + aNc 



AT 

~t7 



AT/Tc < 1, Tc < T. 



(15) 



Therefore, a small increase of the the radius of the ring 5*^ leads to a small increase of the contact term reflected by 
eq. (15). The contact term obviously receives a power like correction ~ e when the geometry is slightly changed ^ e, 

[XYM{T^AT)-XYMiT)] , At AT 

~ e, where e = ~ — — , 

T T ' 



At < T. 



(16) 



Xym{T) ' T 

The correction (16) obviously supports our general argument presented above on power like corrections ~ when 
manifold slightly changes as the variation of the manifold is explicitly translated in changes of the temperature in 
conventional Euclidean formulation. 



* We use term "axion" which is a jargon here. There is no real new dynamical degree of freedom in 4d space such as axion, see reviews [41] 
about the physical 4d axion which is as real degree of freedom and the dark matter candidate. However, the 9 dependent portion of the 
vacuum energy in holographic picture is determined by the dynamics of the axion field living in multidimensional space. Precisely the 
dynamics of the axion field determines the magnitude of the integral (11) which corresponds to the physically observable topological 
susceptibility xym- This relation of the axion field from multidimensional space with XYM justifies our terminology. 
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Now, we want to study a similar variations of the contact term xym {T — AT) with slight change of radius of the 
ring S\ but in confined phase with T ^ T^. This corresponds to a strong coupling computations when no analytical 
formula similar to (14) exists in this regime. From the holographic prospective, however, the the result must be very 
similar to the previously discussed case (16), as the contact term (11) and corresponding vacuum energy, receives 
the power like corrections ~ 'Ylin>\ CpC^j not exponentially small corrections ^ exp(— 1/e) as all changes are governed 
by the massless axion field fieldliving in the bulk of multidimensional space. This argument is consistent with the 
direct lattice computation [44] which indeed finds a power like corrections R~p with p ~ 2.34 for smallest lattice 
size R — 12 and p ~ 1.33 for the largest available lattice size i? = 24 in lattice units. Such a behaviour should be 
contrasted with naive expectation exp(— Aqcd-R) based on conventional dispersion relations without an accounting 
for the non-dispersive term. 

We should emphasize that our general argument docs not specify the power p for the Casimir-like correction ^ e^. 
In general, one could expect any integer number p > 1 as QFT based computations would suggest. Indeed, if 2d QED 
model is defined on a generic torus determined by Teichmiiller complex parameter r = ti + itq the correction behaves 
linearly in inverse size, ^ R^^. The suppression becomes much stronger ^ R^^ (but still, not exponential) in special 
symmetric case when metric t/^i/ ^ diag(l,l) [45]. Also, computations of the vacuum energy in 2d CP^~^ model 
defined on a finite interval of length R with Dirishlet boundary conditions at large also exhibits the correction 
^ R~^ [46]. A similar effect occurs in "deformed QCD" model discussed in section II A when the theory is formulated 
on X R^ and infinitely large space R'^ is replaced by a sphere with radius R. The correction behaves in this 
case as R~^ [47]. One should expect R~^ correction for a symmetric infrared regularization when infinitely large 
space i?^ is replaced by a sphere S'^ [47]. We suspect that a similar variation of power p would also emerge in 
holographic description, however the corresponding study is beyond the scope of the present work. The experience 
with QFT computations suggest that a generic asymmetric background leads to a minimal p = 1 suppression, while 
more symmetric cases lead to a higher power of suppression p > 2. To conclude: all these simple calculable models 
with a gap exhibit the power like corrections, in contrast with naive expectation that there should exponentially weak 
^ exp(— i?) sensitivity to a size of the system R. 

Few more comments: 

We want to avoid any confusion with the terminology. In what follows, we shall use the term "topological Casimir 
effect" to discriminate it from conventional Casimir effect when power like behaviour is due to the real massless 
propagating degrees of freedom living in 4 dimensional space, such as ESzM photons, in huge contrast with our case 
when there are no any massless asymptotic states in QCD. As wc mentioned previously, this "topological Casimir 
effect" occurs in gauge theories when nontrivial "degenerate" topological sectors are present in the system. For the 
present work it is important that the same "topological Casimir effect" in 4d can be thought from the dual holographic 
description as conventional Casimir effect in 5 dimensional space as a result of massless axion field living in the bulk 
of multidimensional space. 

The "topological Casimir effect" obviously is very unnatural and unexpected effect which begs for a simple intuitive 
explanation. A formal explanation was given in the text and was formulated in terms of the topologically protected 
pole at = 0, see eq. (13) within QFT framework. A similar formal explanation in dual holographic description is 
given in terms of the massless axion field living in the bulk of multidimensional space. As this effect plays a crucial 
role in the application considered in section V, we want to present here few other systems where a similar phenomena 
occurs, and where it has precisely the same nature. Furthermore, in these systems a similar problem can be exactly 
solved (in drastic contrast with strongly coupled 4d QCD). Most importantly, analogous effects in these systems have 
been experimentally observed. 

Our first example is from quantum mechanics (QM) and it is the well known Aharonov -Casher effect as formulated 
in [48]. The relevant part of this work can be stated as follows. If one inserts an external charge into superconductor 
when the electric field is exponentially suppressed ~ exp(— r/A) with A being the penetration depth, a neutral magnetic 
fluxon will be still sensitive to an inserted external charge at arbitrary large distance. The effect is pure topological 
and non-local in nature. The crucial element why this phenomenon occurs in spite of the fact that the system is 
gapped is very similar to what we discussed in the present work. First of all, it is the presence of different topological 
states Un (number of Cooper pairs) in the system and "tunnelling" between them (non-vanishing matrix elements 
between ti„ and u„+i states) as described in [48]. Those states are analogous to the topological sectors \n) in our work. 
As a result of the "tunnelling", an appropriate ground state U{6) must be constructed as discussed in [48], analogous 
to the \6) vacuum construction in gauge theories. This state U{9) is an cigcnstatc of the so-called "modular operator" 
which commutes with the hamiltonian. In our work an analogous role plays the large gauge transformation operator 
T such that T\0} = exp{—i9)\9) . An explicit construction of the operator T is known: it is non-local operator similar 
to non-local "modular operator" from ref. [48], see Appendix A for some technical details. The crucial element of 
the construction of ref. [48] is that the induced charges in presence of the gap can not screen the "modular charge" 
as a result of commuting the "modular operator" with hamiltonian. This eventually leads to a non- vanishing effect, 
i.e. collecting a non-zero phase at arbitrary large distance. Wc do not have a luxury to resolve the problem using 
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the hamiltonian description in strongly coupled four dimensional QCD. However, one can argue that the role of 
"modular operator" is played by large gauge transformation operator T which also commutes with the hamiltonian 
[7", H] = 0, such that our system must be transparent to topologically nontrivial pure gauge configurations, similar 
to transparency of the superconductor to the "modular electric field" from ref. [48] . Such a behaviour of our system 
can be thought as a non-local topological effect similar to the non-local Aharonov -Casher effect as formulated in [48], 
see also few comments in Appendix A on this similarity. The last word whether this analogy can be extended to the 
strongly coupled four dimensional QCD remains, of course, the prerogative of the direct lattice computations similar 
to recent studies [44]. 

Our second example is description of a superconductor using the so-called topological "BF" action as presented in 
[49]. This QFT description deals exactly with the same question on how does a moving vortex detect a stationary 
charge, given that the electric field is exponentially screened. Paper [49] explicitly shows how pure gauge (but 
topologically nontrivial) field may lead to long range effects. In many respects this description is very similar to 
our description in terms of unphysical ghost degrees of freedom when the topological features of the system are 
represented by auxiliary fields which are responsible for the dynamics of the degenerate ground state. Our final 
example is description of topological insulators in terms of topological "BF" action as presented in [50] . In this case, 
again, some pure gauge, but topologically nontrivial auxiliary fields may penetrate into gapped insulator to produce 
an effective massless mode on the surface of a sample. Such a behaviour of the system can be thought as a non-local 
topological effect similar to the previously considered example [48]. 

We conclude this section with the following comment. Our interpretation of dynamics of the system, when it is 
described in terms of the different topological sectors \n) and large gauge transformation operator T commuting 
with the hamiltonian [T, i?] = 0, is obviously a gauge variant interpretation, see Appendix A with some technical 
details. As we already mentioned previously, there is no any physical degeneracy in QCD, in contrast with condensed 
matter examples when distinct physical degenerate quantum states are present in the system. Nevertheless, such 
gauge-dependent interpretation helps us to understand a number of very nontrivial features of the contact term which 
are known to be present in lattice Monte Carlo simulations. 

If, instead, we consider a physical Coulomb gauge when all unphysical degrees of freedom are removed from the 
system, the corresponding physics related to nontrivial topological structure of the gauge fields does not go away. 
Instead, it will reappear in terms of the so-called Gribov copies leading to a strong infrared (IR) singularity. Precisely 
this IR behaviour due to the Gribov copies was the crucial element in the approach advocated in ref. [44] where power 
-like correction R^^ has been observed in lattice simulations. Similarly, the construction of the ground state in exactly 
solvable 2d QED in physical Coulomb gauge, as originally discussed in [20], is characterized by long range forces. This 
long range force prevents distant regions from acting independently. This is essentially the same manifestation of 
long range forces in physical Coulomb gauge in the absence of unphysical ghost. We opted to interpret the results 
in gauge-dependent framework in terms of different topological sectors \n) and large gauge transformation operator 
T where one can use analogy with condensed matter systems. Alternative option is lattice numerical simulations 
where entire notion of vacuum winding states \n) and large gauge transformation operator T do not even exist as 
there is a unique ground state. However, irrespectively to the interpretation, the last word in strongly coupled QCD 
is expected from the direct lattice simulations which eventually should confirm or rule out our arguments suggesting 
that the sensitivity to arbitrary large distances should be power like rather than exponential like exp(— Aqcd-R) 
as a consequence of nontrivial topological properties of strongly coupled QCD. 

The most important result of this section which plays a crucial role in the application considered in next section can 
be formulated as follows. The contact term (and corresponding dependent portion of the vacuum energy density) 
may receive power like corrections ^ in QCD despite the presence of a mass gap in the system. From 4d 

viewpoint this is the "topological Casimir effect" when no massless degrees of freedom are present in the system, and 
effect within QFT framework can be explained as a result of degenerate topological sectors in the theory represented 
by topologically protected unphysical pole at k"^ = 0, see eq. (13). The same effect from 5d holographic description 
can be explained in terms of the massless axion field living in the bulk of multidimensional space. From 5d viewpoint 
this is a "conventional Casimir effect" due to a physical (in multidimensional space) massless axion. 



V. CONTACT INTERACTION IN QCD AND PROFOUND CONSEQUENCES FOR EXPANDING 

UNIVERSE 

This portion of the paper which is the direct manifestation of the "topological Casimir effect" discussed in previous 
section is much more speculative in nature than the previous sections. The corresponding speculation on profound 
consequences for cosmology of this effect have been already formulated in author's previous works. However, there are 
two new elements on application of the "topological Casimir effect" to cosmology which were not discussed previously, 
and will be elaborated here. First of all, however, we need to make a short overview of this proposal. 



12 



The 9— dependent portion of the energy which is not related to any physical propagating degrees of freedom, 
is well established effect. It has been tested a numerous number of times in the lattice simulations. How does 
this energy change when the background slightly varies? The main motivation for this question is as follows. We 
adopt the paradigm that the relevant definition of the energy which enters the Einstein equations is the difference 
AE = {E — i^Mink), similar to the well known Casimir effect. Such a definition of the "physical energy" in fact 
is used in description of the horizon's thermodynamics [51, 52] as well as in a course of computations of different 
Green's function in a curved background [53]. In the present context such a definition AE = {E — -EMink) for the 
vacuum energy for the first time was advocated in 1967 by Zcldovich [54], sec also [16] with a large number of related 
references. 

As we mentioned previously, a naive expectation suggests that AE ~ exp(— AqcDZ-ff) ^ exp(— 10*-^) as QCD has 
a mass- gap ~ Aqcd, and therefore, AE must not be sensitive to size of our universe. In this estimate we do not 
distinguish the size of the visible universe ~ from size i? of a compact manifold we used in our analysis in previous 
section. A crucial distinct feature which is important for our discussions is the presence of dimensional parameter 
R ^ H^^ in a system which discriminates it from infinitely large Minkowski space-time. 

This naive argument, however, may fail as a result of nontrivial topological properties of QCD as we discussed 
above. It may lead to a power like scaling AE ^ rather than exponential like AE ^ cxp(— AqcDZ-ff)- From 
holographic viewpoint such power like scaling AE ^ follows from the fact that the contact term (or, what is 
the same, the 9 dependent portion of the energy) is determined by massless axion field. It is naturally to expect 
that massless axion field produces power like corrections in holographic description as discussed in section IV B such 
that surface term (11) which determines the magnitude of the contact term receives a power like corrections when the 
background is slight modified. From holographic perspective, one can view this power like correction as a conventional 
Casimir effect in multidimensional space. However, in our four dimensional space this correction should be interpreted 
as the "topological Casimir effect" as there are no any physical massless asymptotic states in gapped QCD. 

The power like scaling AE ^ with p = 1 (as some explicit QFT computations suggest, see section IV B) may 
have some profound consequences for evolution of our universe. If true, the difference between two metrics (expanding 
universe and Minkowski space-time) would lead to an estimate 

AE ^ i/A^jcD - {lO-^cV)\ (17) 

which is amazingly close to the observed DE value today. Furthermore, the power like scaling (17) with p = 1 would 
imply that our universe approaches the de-Sitter state with constant expansion rate H^c — GAqqj^ in asymptotic 
future [16]. Such a behaviour AE ^ H + 0{H)^ was postulated in [24]. This proposal has received some theoretical 
QFT based support in [16, 21, 47] where it was argued that power like scaling indeed may emerge as a result of non- 
dispersive nature of contact term (3), in contrast with conventional dispersive term (2) related to physical propagating 
degrees of freedom. This correction does not violate unitarity, causality and other important QFT properties. 

We arrive to the same conclusion in the present work using the dual holographic arguments. In the dual picture the 
DE in this framework is due to the Casimir effect in multidimensional space as a result of the dynamics of the massless 
axion field in expanding universe. From four dimensional view point this effect should be interpreted as a "topological 
Casimir effect" when no massless asymptotic states are present in the system. A comprehensive phenomenological 
analysis of this model has been recently performed in [55] where comparison with current observational data including 
Snia, BAO, CMB, BBN has been presented, see also [56-61] with related discussions. The conclusion was that the 
model (17) is consistent with all presently available data, and we refer the reader to these papers on analysis of the 
observational data. 



CONCLUSIONS 

The main result of this work can be formulated as follows. We studied a number of different ingredients related to 
9 dependence, the non-dispersive contribution in topological susceptibility with the "wrong sign" , topological sectors 
in gauge theories, and many related subjects. We argued that the corresponding physics in holographic description 
should be interpreted as a result of tunnelling events governed by the dynamics of the D2 branes. Wc quoted a number 
of theoretical as well as numerical results supporting this interpretation. Essentially, we interpreted the holographic 
formula (11) in terms of a singular behaviour of the contact term (7). The diverging nature of the core and vanishing 
width (in the continuum limit) of this non-dispersive contribution to x is supported by recent lattice studies. After 
all, it is not really a huge surprise that a vanishing width of the contact term (confirmed by the lattice studies) is 
saturated by the objects (also observed on the lattices) which have vanishing sizes and which lead to a superluminal 
behaviour (when interpreted in Minkowski space-time), see original lattice results in [26-34]. 

However, we believe that the most important result of our work is not just a new interpretation of old results [3-5] 
related to the resolution of the U{1)a problem in QCD. Rather, the holographic description allows us to test the 
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sensitivity of the gauge theory with non-trivial topological features to arbitrary large distances. A naive expectation 
based on dispersion relations dictates that a sensitivity to very large distances must be exponentially suppressed when 
the mass gap is present in the system. However, we argued that along with conventional dispersive contribution there 
exists a non-dispersive contribution, not related to any physical propagating degrees of freedom. This non-dispersive 
(contact) term with the "wrong sign" emerges in QFT -based framework as a presence of topologically nontrivial 
sectors and tunnelling events between them. Technically, it is formulated in terms of unphysical massless ghost field 
(13) effectively describing the dynamics of these tunnelling transitions. In dual holographic description the same 
dynamics is described by the massless axion field living in the bulk of multidimensional space. From 4d viewpoint, 
the variation of this contact term with variation of the background leads to a power like "topological Casimir effect" . 
The same effect, but viewed from 5 dimensional perspective can be thought as conventional Casimir effect which is a 
consequence of dynamics of the axion massless field. 

From 4d viewpoint the transparency of the QCD vacuum to topologically nontrivial gauge configurations is similar 
to transparency of a superconductor to the "modular electric field" in the Aharonov -Casher effect as discussed at 
the end of section IV B, see Apendix A with some technical details on this similarity. In both cases an exponential 
suppression of the effects is avoided as a result of topological properties of gauge theories, and can be interpreted as 
a nonlocal effect. We suspect that there should be a holographic description of the Aharonov -Casher effect when its 
"modular electric field" is represented by the axion -like massless field living in the bulk of multidimensional space. 
We leave this subject for the future studies. 

The "topological Casimir effect" in QCD, if confirmed by future analytical and numerical studies, may have profound 
consequences for understanding of the expanding universe we live in. Finally, what is perhaps more remarkable is the 
fact that some elements of this framework can be, in principle, experimentally tested in heavy ion collisions, including 
the Casimir like L~p scaling behaviour. The parameter L in heavy ion collisions describes a size of the region where 
the QCD vacuum is disturbed as a result of collision. Dependence on L, in principle, can be experimentally studied 
as L is related to the physical size of the colliding ions, see [25, 62] for the details. 
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Appendix A: The large gauge transformation operator T. 

The goal of this Appendix is to give a short overview of properties of the T operator which was essential element in 
our discussions of the physical interpretation of the contact term. We also want to present some additional arguments 
demonstrating striking similarity between T operator and "modular operator" in the Aharonov -Casher effect as 
mentioned at the end of section IV B . 

The basic properties of the T operator were analyzed in late seventies and we follow classical paper [63] in description 
of some features of this operator in the context of the present work. The starting point is the construction of the 
unitary operator T, 

r = e'«^ , Qx = J d'xA'l (V,A'^ + fabcAtX"") , (Al) 

where we literally use notations of ref. [63]. An arbitrary function A''(x) in this formula is a c- function which is 
a parameter of the gauge transformations. In particular, if X'^^x) vanishes at spatial infinity one can represent the 
operator Qx in form 

= - y d^xX^ix) (V,i^ + /aOcA^i?) = - y d^xX^{x)C''{A), (A2) 

which obviously annihilates all physical states as the term in parentheses is just Gauss's law C{A)\7iphys) ~ 0. 
However, if X°'{x) describes the large gauge transformation, the operator T becomes nontrivial operator and describes 
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the transition between topologically different sectors of the theory: 

T\n)^\n + l). (A3) 
In this case one should construct the so-called \9) vacuum state which is an eigenstate of the T operator: 

|0)=^e™V), m^e-'^y). (A4) 

n 

It is important to emphasize that while operator T formally constructed as an operator of gauge transformations, 
this operator does change the state (A3) as a result of global effect. Therefore, one should treat T as "improper" 
gauge transformation (the "large gauge transformation"). Still, T commutes with the hamiltonian [T, if] = 0. As 
we mentioned in the text, there is no any physical degeneracy in QCD, in contrast with condensed matter examples 
when distinct physical degenerate quantum states are present in the system. Nevertheless, the gauge-dependent 
interpretation in terms of the tunnelling between |n) and l^^-l- 1) states is quite useful to develop an intuitive picture of 
relevant physics, which otherwise hidden in the numerical lattice Monte Carlo simulations with such puzzling elements 
as the "wrong sign" contributions with vanishing width as discussed in section IV. 

While this paper is mainly devoted to 4d QCD it might be instructive to present the large gauge transformation 
operator T^d for 2d QED as well. In this case, as it is known, the nontrivial topological structure emerges as a result 
of nontrivial mapping 7ri([/(l)) = Z which replaces the classification of the topological sectors in QCD which was 
based on ■k^{SU{N)) = Z mapping. In 2d case a similarity between the "modular electric field" in the Aharonov 
-Cashcr effect and operator 72d becomes much more apparent. Indeed, the large gauge transformation operator 72d 
for 2d QED can be represented as follows 

T2d = e-'ir^^'-) . ei f-- '^^K^)ciA) ^^^^ 

where g is the coupling constant in 2d QED, and E ~ —A^ is electric field. Formula (A5) was obtained by integrating 
by parts the two-dimensional analog of eq. (Al) and taking into account that the large gauge transformation in this 
case is characterized by X{x = oo) — X{x = —oo) = l-n. One can construct \Q) state similar to (A4) with result that 
the expectation value of large gauge transformation operator Tid is given by 

m2d\Q)=e-'\ e=—E,{oo). (A6) 
9 

The fact that in 2d QED the electric field at infinity plays the role of the 6 parameter of course is well known classical 
result [64]. A similarity of the expectation value (A6) with the Aharonov -Casher effect as presented in ref. [48] is 
striking as 9 parameter in this case is almost identical to the "modular electric field" from ref. [48], with the "only" 
difference that \n) states in 2d QED are not physically distinct degenerate states as we already emphasized, in huge 
contrast with physical degeneracy in condensed matter systems discussed at the end of section IV B. Eventually, this 
difference leads to the "wrong sign" in topological susceptibility as we argued in this work. The last word whether 
this analogy can be extended to the strongly coupled four dimensional QCD remains, of course, the prerogative of 
the direct lattice computations similar to recent studies [44]. 
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